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Abstract 

A possibility of describing two-level atom states in terms of pos- 
itive probability distributions (analog to the symplectic tomography 
scheme) is considered. As a result the basis of the irreducible represen- 
tation of a rotation group can be realized by a family of the probabil- 
ity distributions of the spin projection parametrized by points on the 
sphere. Furthermore the tomography of rotational states of molecules 
and nuclei which can be described by the model of a symmetric top is 
discussed. 



Recently JT|, |2| the symplectic tomography scheme was suggested to obtain 
the Wigner function by measuring the probability distribution function (the 
marginal distribution) for a shifted, rotated and squeezed quadrature, which 
depends on extra parameters. The invertable map of the Wigner function 
of a quantum state onto the positive marginal distribution for the continous 
observable (coordinate) was used to give the formulation of quantum dynam- 
ics as the classical statistical process H [|, |5|. From this point of view, the 
Moyal approach || to the quantum evolution as to a statistical process was 
improved in the sense that instead of the Moyal quasidistribution function 
(the Wigner function) was introduced the positive probability distribution 
of the measurable variables describing an arbitrary quantum state and its 
evolution. A spin states reconstruction procedure similar to the symplectic 
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tomography scheme is investigated in J7L M. The authors of these papers 
showed that in the framework of the symplectic tomography approach it is 
possible to describe equivalently the spin state in terms of the classical dis- 
tribution function of a discrete variable instead of the wave function or the 
density matrix. Taking into account the results of these papers in this article 
it will be investigated the tomography scheme for two-level atom states and 
rotational states of molecules and nuclei. 

The goal of this work is to make a review of construction of the explicit 
formula for the invertable map connecting the quantum state of a two-level 
atom described by the 1/2-spin density matrix p^ 1 / 2 ) with the probability 
distribution function w (±1/2, 9, ip) of the 1/2-spin projection on the quan- 
tization axis, where the angles < # < 7r, < ip < 2tt determine the 
coordinates of the vector normal to the surface of the sphere of the unity ra- 
dius. On the other hand, quantum characteristics such as rotational energy 
levels of molecules and nuclei are connected with the tomographic descrip- 
tion of the tops' behavior. Therefore in this work it will be also discussed 
the classical-like description of the quantum states of a top in terms of the 
distribution functions. For simplicity of presentation, we will consider the 
example of a symmetric top. These results also are discussed in ||, |H], [TT| . 



1 Tomography of two-level atom states 



It is a common knowledge (see, for example |12|]) that the wave function ip (m) 
of the 1/2-spin particle, which can describe the state of a two-level atom and 
consists of two components if) (1/2) and if) (— 1/2), can be represented in the 
form of a spinor 

*-(#$>)-(;)■ 

In the case of the pure particle's state, the Hermitian density matrix has the 
form of a 2 x 2 matrix 

,»/2> = ^ = ( M 2 f i2 ) i (2) 
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where the diagonal matrix elements satisfy the normalization condition | a 
+ | b | 2 = 1. Any rotation in the three-dimensional space, determined by 
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< < 2tc, is associated with a 2 x 2 unitary matrix 
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Under rotation, the Hermitian density matrix 

p (l/2) is 

transformed as 

p (l/2) ^ p (l/2) (m) = up (l/2) utj ( 4 ) 

where the 2x2 unitary matrix m is a representation of the angle that de- 
scribes the matrix elements of the matrix u. For the diagonal elements of 
the Hermitian density matrix p^ 1 / 2 ) (u) we get 



pZ /2) (u) 



1/2 

E 



1/2 

E 

l/2m=-l/2 



A ( s 1/2) HP^ 2) A ( i /2) *W 



(5) 



Here -D^ 2 ^ is the Wigner D- function and i = —1/2, 1/2. The diagonal 
elements of the density matrix of the quantum state take nonnegative values, 
and their sum is equal to unity. The physical meaning of these elements is 
that they are a probability to measure the value of the spin projection on 
the fixed axis in a space. Then, we introduce the notation 



p!* /2) («) 



w 1,11 



(6) 



where the function w(i,u) is the marginal distribution, i.e., the probability 
of finding the spin projection i on the fixed axis in the space for the 1/2- 
spin particle. From Eq. @, we see that the marginal distribution also 
depends on the Euler angles <p, 9, ip as parameters. This distribution function 
is normalized for all values of the Euler angles. From the structure of Eq. 
(H), it follows that the marginal distribution w(i,u) depends only on two 
Euler angles, and does not depend on the angle ip of the rotation. In fact, 
according to Eq. (|), one can obtain the marginal distribution w (i, u) of the 
1/2-spin particle in the following form 



w(l/2,u) 



2 9 , l2 sin# iu} sin^ 
cos 2 - a 2 H e lv ab H 



e~ lLp ha* + sin' 
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(7) 
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and 

o 9 , ,o sin 6 .• . . sin 6 •„, . o 6 . , l0 , . 

w -1/2, it) = sin 2 - a 2 e^afo* e^fta* + cos 2 - 6 2 . 8 

2 1 1 2 2 2 1 1 w 

Hence, the family of the probability distribution functions of the 1/2-spin 
projection is parametrized by the point's coordinates 0, if on the sphere of 
unity radius. This parametrization coincides with the physical meaning of 
the marginal distribution in the sense that the distribution function w (i, u) 
is the probability to observe the spin projection i if we measure this spin 
projection on the quantization axis which is parallel to the vector normal to 
the surface of the sphere of the unity radius in the point with the coordinates 
9 and <p. If we know the positive, normalized marginal distribution w (i, u), 
then, as it was shown in [0, |S|, the matrix elements p^ m , can be calculated 
with the help of the measurable marginal distribution w (i,u) of the particle 
with an arbitrary spin j and the values of indices i = —j, + by 
means of the relation 

(-i) m 'p!L = E E (2A: + i) 2 E (-iy 

k=0 l=-k 





w (i, u) D Ql (u 

where m, m! = —j, + and the integration leads over the rotation 

angles f,9,ip 

JdQ = j df j d^j J sm6d6. (10) 



In the case of the 1/2-spin (i.e., j = 1/2), Eq. (H) reads 

p%$ = (-ir'EE(2Hi) 2 E (-1/ 

k=0l=-k i=-l/2 

w(i,u)D k a {u)—A \ 2. k )( \ 3 

87r 2 \ % -i J \m -ml I 

where m, ml = —1/2, 1/2. Calculating all non-zero integrals with D- 
functions and the marginal distributions and taking into account the ex- 
pressions of 3j-symbols (see, for example, [0), we can rewrite Eq. (|11|) as 
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follows 

+ * ( 5 5 j)(|„|'-|M') 
y/Q V m — m J V 1 1 1 1 / 

W3( 5 5 1 )ab*-iVS{ 5 2 j W, (12) 
V m — m — 1 / \ m — m 1 / 

where m, m' = —1/2, 1/2, and i is the imaginary unity. The substitution 
of m, m! = —1/2, 1/2 in Eq. (|T2p leads to the following result: 

Pmm> =1 a m = m'= 1/2, p^J =| 6 | 2 , m = m'= -1/2, (13) 

and in all the other cases 

Pw m = -m' = l/2, p^ = ba% - m = w! = 1/2. (14) 

Then, we checked by the direct calculation that the right-hand side of Eq. 
([TT|) is equal to the density matrix fl2|). One can conclude, that given a 
measurable marginal distribution of a particle, whose state is described in 
terms of spinor, one can reconstruct the state density matrix by means of 
Eq. 



2 Examples of the marginal distribution for 
the 1/2 and 1 spin states 

The classical-like description of quantum mechanics in terms of the positive, 
normalized marginal distribution can be easily understood in the case of the 
1/2-spin particle. It is a common knowledge |12"| |, that the spin operator of 
the 1/2-spin particle has the form 

s=h, (15) 

where the Pauli matrixes are 

d - = f i n J > °v = ( _,• ~n ) > ^ = ( n _1 J • ( 16 ) 
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The wave function of the 1/2-spin particle can be represented in the form 
of spinor ([!]). When the direction of the particle's spin coincided with the 
positive direction of the x axis in the coordinate system (x, y, z), the wave 
function of spinor has the form 

e'^(l). (it) 

And in the case of a contrary direction of the spin with respect to the x axis, 
the wave function reads 

= 75 ( -i ) " (18) 

Analogously for the ^-projection, when the direction of the particle's spin is 
+ 1/2 or —1/2 on the y and z axes, we have 

*s +, =4(lV us) 



y/2 V i 
1 

V2 



and for the z-projection 



W = ( I j , (21) 

= ( ; ) • (22) 

For all values of the 1/2-spin directions on the x, y and z axes the density 
matrices have the following forms 

& = \ !)■ & = \{-i (23) 

<4 + H0 t). *>-§(i 0- (24) 

and 
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Applying Eqs. (0) and (El), we obtain the marginal distributions for the 
different values of the 1/2-spin directions on the x, y and z axes 

wl +) (+1/2, u) = i (1 + sin cos y?) , (-1/2, u) = ]-(!- sin cosy?) 

(26) 

for the positive x-projection, 

w i ~ ) (+l/2,u) = - (1 - sin 6 cosy?) , w { ~ ] (-1/2, u) = - (1 + sin 9 cosy?) 

(27) 

for the negative x-projection, 



+1/2, u) = - (1 + sin6»siny?) , w/ +) (-l/ 2 ,«) = - (1 - sin6»siny?) 



2 

for the positive y-projection 



(28) 



for the negative y-projection 

- ens 2 

2 



w y ~ ] (+1/2, u) = - (1 - sinflsiny?) , w y ~ ] (-1/2, u) = - (1 + sin6>siny?) 

(29) 



(+1/2, it) = cos 2 -, «4 +) (-1/2,m) =sin 2 - (30) 



for the positive ^-projection and 

w^{+l/2,u) = sin 2 ^, u;H(-l/2,«) = cos 2 ^ (31) 

for the negative z-projection. In Figs. 1-2, we plot the marginal distributions 
of the 1/2-spin particle (±1/2, 9, y?) the spin direction of which coincides 
with the positive direction of the x axis and the spin projection onto a fixed 
axis in space is equal to ±1/2, as a function of the angles 9 and y?. 

However, there exist the mixed states of the 1/2-spin particle (see, for 
example, WW )- These states are described only by the density matrix 

*-».)- (32) 

In Eq. (j32D , s± are determined by the relation s± = s x ± is y , and s x , s y and 
s z are the mean values of the 1/2-spin projections on the x, y and z axes, 
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respectively. These parameters satisfy the condition s 2 x + s 2 + s 2 z < 1 / A. A 
parameter \i = Tr (p^J is named "purity state degree" and depends on the 
mean values of the spin projection p, = 1/2 + 2 (s 2 + s 2 + sf\ . As in the case 
of a pure state, using Eq. (|5|) for the matrix elements of the density matrix 
of a mixed state, we calculate the marginal distribution of the mixed state 
of the 1/2-spin particle, 

w m (- , u^j = - + s z cos 9 + s x sin 9 cos ip + s y sin 9 sin ip (33) 

and 

u> m ^— - , M^j = - — s z cos — s x sin 6* cos (p — s y sin sin (p. (34) 

It is easy to check that the examples of the pure-state marginal distribution 
(|26D-(j31~|) are obtained from Eqs. (|33|), (|34l) by the proper choice of the 



parameters s x , ~s y and s z . 

In this section, one can also consider another one important example of 
the 1-spin particle states, the direction of which coincides with the positive 
direction of the z axis. In this particular case, the wave function has the 
form 

W = ( j , (35) 

and the density matrix, consequently, takes the form 

/ 1 \ 
= 0. (36) 
V / 

The marginal distribution of the spin state can be calculated by the following 
equation (see, for example, |fP2"| ) 

pi 1} (u) = D\ t («) P V D * (u) , (37) 

where u is determined by Eq. (p~0|) . As a result, we obtain the marginal 

distribution w (i, ip, 9, ip) = (u), Z)j=-i w {h 0, ^, v 9 ) = 1? i- e -> three proba- 
bilities for the 1-spin projection onto the z axes (—1, 0, and +1), which also 
depend on the rotation angles of the reference frame <p, 9, and ip 

, r\ \ (1 + cos 9) 2 , , 

w(l,9) = ± (38) 
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w(0,6) 



cos 



(39) 



COS ( 



W I 



(40) 



A dependence on the 0, (p angles drops out. Substituting this marginal 
distribution into formula (^) by means of which matrix elements of the density 
matrix p^ m i can be reconstructed by the measurable marginal distribution 
w(i,u), i = —1,0,1 and executing the calculations, we could verify the 
correction of this formula in the particular case of the unity spin. In the 
case j = 1, this formula consists of twenty seven items. Let us consider all 
values of the Wigner D-functions for the various values of their subscripts 
and superscripts. Because the marginal distribution w(i,9), i = —1,0,1 
takes only positive values, and because some D-functions contain a complex- 
valued exponential coefficient, integration of the product of the marginal 
distribution and a D-function over the 0, 8, ip rotation variables yields zero. 
Taking this fact into account, Eq. (pi) can be simplified, and we arrive at 



o {1) 

rmm' 



("l) r 



E fdn (-iyw(i,e)D° 00 



+9 E Jdn(-l) l w(t,6)Dl 



oo 



+25 E Idtt{-iyw(i,6)D 



00 



1 1 0W1 1 
i —i J I m —vnl 
1 ' 1 1 \ / 1 1 1 \ 
% —i J I m —m! J 
1 1 2 \f 1 1 2 \ 
i —i \ m —m' / 



(41) 



Calculating all nonzero integrals and taking into account the expressions of 



the 3j-symbols from |12j we can rewrite Eq. ([II]) into the form 

-m'+l 



o {1) 

rmm' 



1 





(42) 



For m m! (m, m' = —1, 0, 1), 3j-symbols are equal to zero (see, for 



example, 
obtain 



Evaluating all the remaining 3j-symbols for m 



m , we 



o (1) 

rmm 1 




-1, 0, 1. 



(43) 
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This means that we demonstrated how the initial density matrix fl36|) can 
be recovered, if we know the marginal distribution of the quantum state 



h 2 c ft 2 / 1 1 



H = w: J+ T{r + rJ J f < 44 > 



3 Tomography of rotational states of 
molecules and nuclei 

Here we will consider the tomography scheme of rotational states of molecules 
and nuclei using, for simplicity, the model of a symmetric top. It is a common 
knowledge (see, for example, |l!|) that the Hamiltonian of a symmetric top 
has the form 

t-.-' 

J 2 

where I a, Ic are the principal inertia momenta of a top (two of the mo- 
menta of a symmetric top coincide with each other), and J is the angular- 
momentum operator. Stationary states of a symmetric top are characterized 
by three quantum numbers: the orbital momentum j and its projections onto 
the top axis (J^ — k,k — — j, — j + 1 . . . j) and onto a fixed z axis in space 
(J z = M, M = —j, — j + 1 . . . j). The symmetric top energy is independent 
of the last quantum number M. Let us consider the stationary states of 
the symmetric top with given energy. To do this, we take a subspace of the 
(2j + l) 2 dimensions with the fixed j in the state space. The wave function 
of the stationary state of a symmetric top can be represented in the form 

Vfii = (jMk | ^> = ^$D{ M , (45) 

where D J kM is the Wigner D-function, and ipjj) is the wave function in the 
reference frame, which comoving with the physical system (the top). For a 
pure state, the density matrix of the symmetric top is expressed in terms of 
the wave functions as follows 

PMkM'k' = i'Mki'M'k'- (46) 

Thus, under two consecutive rotations (these rotations are determined by the 
Euler angles u((j),6,<p) and u' (tp 1 ,8' the Hermitian density matrix p is 
transformed as 

p -> p {j) (u, u') = D {j) («') D {j) (u) pWDW (u) DWt ? (47) 
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where unitary (2j + 1) x (2j + 1) matrixes D (u) and D (u 1 ) represent these 
two rotations through the angles u and u 1 . Since the density matrix under 
consideration depends on four discrete indices, it is necessary to do two con- 
secutive rotations of the reference frame, which comoving with the top to 
obtain the diagonal elements of the transformed density matrix. Taking into 
account this fact, for the diagonal elements of the density matrix p& (u, u') 
can be reduced to the form 

Pgrfxfc («, «0 = I>U M D nn («) P%A («) («') • (48) 

Here -D^ (it) is the above Wigner D-function and ii, i 2 , n,p, s, I = —j, —j + 
l,...,j. In Eq. (|8|), we assume that summation is performed over the 
repeated indices n, p, s, I, but there is no summation over the indices z'i, 
Below, we will assume that the Hermitian nonnegative density matrix p£, s ; 
describes not only the pure state of the top fl45|) but also an arbitrary mixed 

state, i.e. Tr(p^A < 1. 

Let us discuss the problem of reconstructing all matrix elements 
the diagonal elements Pi li2ili2 (u, u') are known. To do this, let us multiply 
both sides of Eq. ( fl8|) by the functions D c ab (u) and D{ e (u') with arbitrary 
upper indices c, /, and with different angular variables u and u' . As a result, 
we have 



pfnniA^)Dl b {u)D{ e (u') = Dl l2 {u')Dl n {u)p^ sl Dll{u)D{: 2 



hi 



®D c ab (u) Di {v!) . (49) 
Integrating over the angular variables u (0, 9, if) and u' (<j)', 6', ip'), 



27T 7T 27T 

1 



J smOdO J d(p, (50) 





2-7T tt 2n 



J dtt = i j d<j>' J sin 9'dO' J dip', (51) 





using the well-known expression of the integral of the three D-functions prod- 



uct IJ] and 

D 3 : /m (u) = (-l) m '- m Dl m ,_ m (u), (52) 
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along with the orthogonality properties of the 3j-symbols, one can express 
p^ipsi in terms of the diagonal elements p^i 2ili2 (u,u ! ). As it was mentioned 
above, the diagonal elements of the density matrix for a quantum state take 
nonnegative values and their sum is equal to unity. Thus, we introduce the 
notation p^] 2il j 2 (u, u') = w (ii, i 2 , u, u'), where the function w (ii, i 2 , u, u') is 
the marginal distribution of a symmetric top, i.e. the probability of finding 
the projection i x of the angular momentum onto a fixed axis in space and 
the projection % 2 of the angular momentum onto the top axis, and also this 
probability parametrically depends on Euler angles cj),9,ijj,(j)',9',i/j'. This 
function is normalized J2i!,i 2 =-j w (*i> ^2, u, u') = 1 for all values of the Eu- 
ler angles. Suppose we know the positive, normalized marginal distribu- 
tion w (ii, %%, u, u'). Thus, in view of Eqs. p5|)-([52]) the matrix elements 



U) 

mil 

tion w (ii, i 2 , u, u') with i\, i 2 = —j, — j + 1, . . . , j by using the relationship 



Pm 1 m 2 m'rn' 2 can ^ e expressed in terms of the measurable marginal distribu- 



2j ki 2j k 2 

(2k 2 + lf 



p { £ imM = E E (2*1 + 1) 2 E E (2* 2 + i 




w (ii,i 2 ,u,u') DqI (u) 



—m\ 



(53) 

Here mi, m 2: w! 1: m' 2 = + One can conclude that given a 

measurable marginal distribution of the symmetric top, one can reconstruct 
the stationary-state density matrix of the symmetric top by means of Eq. 
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Figure 2 




